NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  -vdiatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  fonmilated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  Euiy 
other  person  or  corporation,  or  conveying  any  ri^ts 
or  pemlsslon  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


FOREIGN  TECHNOLOGY 
DIVISION 

AIR  FORCE  SYSTEMS  COMAAAND 


WRIOHT-PAnEIISON  ANt  FO*a  BASE 
OHIO 


FTD-TT-  62-1694/1+2 


UNEDITED  ROUGH  DRAFT  TRANSLATION 

i  . 

I 

I. 


JOURNAL  OF  ENGINEERING  (SELECTED  ARTICLES) 
Eaglish  Pages ;  43 

SOURCE;  Russian  periodical,  Inzhenernyy  Zhurnal, 
•  Vol.  2,  Nr.l,  1962,  pp.  79-86,  163-170, 
175-181 


'  S/258-62-2-1 


THIS  TRANSLATION  IS  A  RENDITION  OP  THE  ORlGL 

HAL  FOREIGN  TEXT  WITHOUT  ANY  ANALYTICAL  OR 
EDITORIAL  COMMENT.  STATEMENTS  OR  THEORIES 

PREPARED  BYi 

- 

ADVOCATED  OR  IMPLIED  ARE  THOSE  OF  THE  SOURCE 

AND  DO  NOT  NECESSARILY  REFLECT  THE  POSITION 

TRANSLATION  SERVICES  BRANCH 

OR  OPINION  OF  THE  FOREIGN  TECHNOLOGY  OL 

PORIICH  TECHNOLOOY  DIVISION 

• 

VISION*  1 

WP>APB,  OHIO. 

FTD-TT-  62-1694/1+2 


Date  29  March  19  63 


S 


TABLE  OF  CONTENTS 


A.  I.  Tolstykh,  Turbulent  Boundary  Layer  with  Pressure.  Gra¬ 
dients  on  a  Porous  Surface . .  1 

V.  A.  Bronshten,  The  Processes  of  Ionization  and  Recomblna-^ 

tlon  In  Hot  Air . .  .  .  ■ . .  ,  '  15 

Yu. I.  Krasil'nikov,  Aerodynamic  Characteristics  of  Delta 

Wing  with  Deflected  Aileron,  Flap,  and  Other  Secondary  Con- : 

trol  Surfaces  In  Supersonic  Flow  .  .  . . .  29 


-  1  - 


PTD-  TT-62- 1694/1+2 


TURBULENT  BOUNDARY  LAYER  WITH  PRESSURE'  GRADIENTS 
ON  A  POROUS  SURFACE 
A.  I.  Tolstykh 
(Moscow) 

Existing  methods  of  calculating  a  turbulent  boundary  layer  with 
the  transport  of  a  substance  through  the  surface  are  restricted  only 
by  the  case  of  zero  pressure  gradients  and  are  based  on  the  use  of  the 
same  semlemplrlcal  relationships  as  in  the  case  of  an  impermeable 
plate  [l]-[4],  [8],  Because  of  the  lack  of  the 'required  experi¬ 

mental  data  the  latter  unavoidably  results  in  the  assumption  that  cer¬ 
tain  quantities  are  Independent  of  the  mass -transfer  parajneters  and  to 
a  somewhat  arbitrary  approach  in  the  selection  of  the  necessary  con¬ 
stants. 

The  purpose  of  the  present  work  is  to  investigate  the  more  gen-, 
eral  case,  l.e.  ,  a  turbulent  layer  on  a  porous  surface  in  the  presence 
of  pressure  gradients,  without  the  introduction  of  any  semiemplrlcal 
relationships  and  with  the  utilization  of  only  that  information  which 
can  be  obtained  from  the  known  prior  history  of  the  flow. 

Initially  an  incompressible  boundary  layer  is  considered  in  the 
absence  of  diffusion,  and  then  a  formal  generalization  is  carried  put 
for  the  case  of  a  heterogeneous  mixture  of  compressible  gas. 

1.  To  derive  the  required  relationships  we  will  proceed  from  the 
case  of  an  idealized  two-layer  model  and  will  assume  a  turbulent  layer 
consisting  of  a  wall  region  in  which  the  equations  of  laminar  motion 

1  _ 
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and  a  turbulent  core  are  valid.  In  accordance  with  this,  one  part  of 
the  flow  may  be  determined  by  assuming  the  required  number  of  boundary 
conditions,  and  with  respect  to  the  other  part  of  the  flow  (because 
the  equations  of  the  turbulent  boundary  layer  are  open)  it  becomes 
necessary  to  make  certain  assumptions. 

Flow  in  the  wall  region  (a  viscous  sublayer)  within  the  scope 
of  the  boundary -layer  equations  may  be  described  by  a  system  of  the 
following  form; 


du 


du 


ax  dy  ’  / 


(1) 


with. boundary  conditions  at  y  =  0 

«=o,  »,=o,  = 

•  Vp  1 

Here  v^  =  v  —  w;  u  and  v  are  velocity  components  In  the  conven¬ 
tional  boundary -layer  coordinates  x  and  w  is  the  normal  velocity  at 
the  surface]  p,  and  v  is  the  stress  at  the  wall,  and  the  density 

and  kinematic  viscosity]  and  are,  respectively,  certain  func¬ 
tions  characterizing  the  flow  and  their  value  at  the  wall.  For  exam¬ 
ple,  in  the  presence  of  a  pressure  gradient  =  -( l/p)  (dx/dp) , 

as  a  gas  is  passed  through  the  surface  ZQ^  "  ^2  ~ 

At  given  values  of  and  System  (l)  fully  determines  the  mo¬ 
tion  near  the  wall  and  in  view  of  the  smallness  of  the  nonlinear  mem¬ 
bers  is  reduced  to  the  following  equation: 

dy* 


with  the  corresponding  boundary  conditions.  Prom  the  concept  of  dimen¬ 
sionality  it  follows  that  the  solution  to  the  above  equation  has , the  ■ 
following  form: 
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where  =  V^(t  /p)  is  the  dynamic  velocity. 

*  w  . 

In  particular,  when  air  Is  forced  through  a  surface  along  which 
there  is.  the  following  pressure  gradient: 


r 


where  y^  =  (yv^/v),  and  a  and  f  are  the  dimensionless  parameters  of 
the  pressure  gradient  and  the  air  injection: 


a* 


i  dp  V 

P  dx  * 


(2) 


The  wall  boundary  layer  in  which  the  flow  Is  described  by  a  "lam¬ 
inar"  solution  Is  naturally  defined  as  a  line  along  which  laminar 
friction  covers  a  definite  and  fully  established  portion  of  the  turbu¬ 
lent  layer;  let  us  say  that  It  Is  defined  as  the  line  along  which 
v^/v  =  const,  where  Is  the  turbulent  (virtual)  viscosity  that  Is 
associated  with  turbulent  friction  by  the  relationship  =  pv^(du/dy). 
This  representation,  accurate  to  the  constant,  defines  the  "laminar 
sublayer"  on  the  Impermeable  piate,  since  the  relationship  v^/v  must  > 
be  a  function  only  of  y^^.  In  the  general  case,  this  relationship  may 
be  a  function  also  of  the  parameters  However,  If  instead  of 

the  local  value  of  total  friction  t  Is  taken  as  one  of  the  determining 
parameters,  the  effect  of  must  make  Itself  felt,  apparently,  only 
by  changing  the  value  of  t.  Indeed,  on  the  one  hand  t  =  pv(1-{t— ) — ,  . 

■  V  /  dl/  .  —• 

while  on  the  other  hand,  as  a  result  of  evaluating  the  various  terms 
in  the  balanced  equation  for  pulsating  energy  [10],  (v^/v) (bu/Zy)^  we 
(e  denotes  those  terms  which  correspond  to  the  averaged  viscous  dissi¬ 
pation  of  pulsation),  only  If  convection  transfer  of  the  pulsating  en- 
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ergy  by  the  lateral  velocities  Is  not  too  great  (if  the  injection  ve¬ 
locity  w  is  not  too  great)  and  y*  is  not  too  small* 

In  view  of  the  fact  that  e  is  not  an  explicit  function  of  and 
is  most  likely  determined  by  t  (by  changing  the  velocity  scale  for 
major  pulsations),  the  ratio  v^/v  is  a  function  only  of  t,  p,  and 
v,  1. e. ,  the  condition  v^/T  =  const  is  equivalent  to  the  relationship 
(y/v)>/(T/p)  =  const.  Since  in  an  idealized  two-layer  model  the  total 
stress  near  the  wall  is  only  viscous,  we  require  that  at  the  boundai*y 
of  the  viscous  sublayer  the  following  relationship  be  satisfied: 


V 


(3) 


where  the  constant  in  the  right-hand  part  is  independent  of  the  x  co- 

@ 

ordinate  and,  consequently,  it  is  Independent  of  the  parameters  of  In 
Jectlon  and  the  gradient.  This  constant  can  be  determined  if  the.  re¬ 
quired  magnitudes  pertaining  to  some  "initial"  section  are  known;  in 
particular,  if  the  boundary  layer  on  the  Impermeable  plate  has  been 
thoroughly  studied  in  the  prior  history  of  the  given  flow,  const  =  a, 
where  a  is  a  universal  constant.  In  view  of  the  universality  of  d  we 
can  maintain  that  the  following  is  valid  regardless  of  the  existence 
of  such  a  prior  history: 

*  litn  I  /  V =  const  =«, 

rZ!  y 

Thus  Condition  (3)  can  be  written  in  the  following  form: 


K*JL  (4) 

where  Y  is  the  value  of  y^  at  the  boundary  of  the  viscous  region,  and 


Let  us  now  examine  the  turbulent  portion  of  the  boundary  layer. 
Restricting  the  classes  of  flows  being  studied  somewhat,  we  will  as- 
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sume  that  the  velocity  profile  In  any  section  can  be  presented  In  the 
"single -parameter"  form 

,  ~  f  [t’  "  ’  f  =  *• , 

where  6  Is  the  thickness  of  the  boundary  layer;  Is  the  velocity  at 
the  outer  limit  of  the  boundary  layer;  r  Is  some  parameter.  In  the 
"Initial"  section  the  velocity  distribution  must,  thus,  be  described 
by  a  function  of  the  form  of  (5)  at  some  Initial  value  of  r.  In  order 
to  clarify  the  nature  of  the  function  we  will  use  the  concepts 

of  dimensionality.  The  only  magnitudes  determining  the  flow  In  the 
vicinity  of  the  boundaries  of  these  regions  will,  evidently,  be  the 
total  momentum  flux  t  at  this  boundary,  density  p,  and  the  distance 
from  these  parameters  we  can  make  only  the  single  combination  having 

'l/Z 

the  dimension  of  the  derivative  du/By,  and  namely  =  const 

y  ..I. 

therefore  at  the  boundary  as  a  consequence  of  Relationship  (3)  and.  the 
condition  of  momentum-flux  continuity  we  find  the  following  equality 


du  .  &v^q 

—  =  const  — 
dy  dy 


(6) 


which  Is  accurate  to  the  higher  derivatives  of  the  velocity  u  and 
yields  the  function 

With  respect  to  the  form  of  the  function  f  let  us  note  that  For¬ 
mula  (5)  must  determine  the  velocity  profile  even  In  the  case  In  which 
the  "Initial"  section  corresponds  to  the  flow  on  the  Impermeable 
plate;  It  Is  evident  that  Assumption  (5)  will  be  satisfied  if  the;  cor¬ 
responding  functional  relationship  will  have  the  same  form  as  In  the 
Initial  section,  l.e. ,  for  example. 


In  x  oi:JL^/JLr 

*(<?<•)  •  *  '  v«/ 


(7) 


where  the  parameter  r,  denoted,  respectively,  by  k  and  m,  assumes  the 
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known  values  k(0)  =  and  m(0)  =  itIq. 

As  shown  by  the  experimental  data,  the  second  of  the  functions 
(7)  describes  an  extremely  extensive  class  of  flows  (in  particular, 
flows  with  . nonpositive  pressure  gradients  and  flows  with  the  movement, 
generally  speaking,  of  a  foreign  gas  through  the  porous  surface  of  the 
plate  [5-7] )•  It  Is  good  that  the  step  distribution  of  velocity  Is 
preserved  to  substantial  Injection,  Intensities  [5],  [7].  Subsequently 
we  will  use  the  step  function  which  results  In  comparatively  simple 
finite  results;  the  exponential  values  will  be  found  from  Relation¬ 
ship  [6],  In  which  the  constant,  as  a  result  of  the  universality  of 
the  Initial  value  of  mQ,  Is  determined  as  a  result  of  the  limit  tran¬ 
sition  for  “*•  0:* 

.  It  (8) 

In  the  derivation  of  this  formula  we  used  the  condition  of  ve¬ 
locity  continuity  which  can  also  be  written  In  the  following,  foiin: 


where  (c^/2)  -  (v^/u^)  and  R  =  (Ug6/v).  Relationships  (4),  (8),  (9), 
and  (2)  make  It  possible  to  calculate  the  local  coefficient  of  fric¬ 
tion  and  the  exponent  m,  which  when  using  the  Integral  momentum  .^rela¬ 
tionship  fully  solves  the  problem. 

.  ,  Let  us  examine  certain  special  cases. 

a)  Let,  for  example,  f  =  0,  and  a  <  0  (the  negative  pressure- 
gradient  on  the  Impermeable  wall).  In  this  case  the  velocity  profile 
has  a  step  shape,  and  Formulas  (4)  and  (8)  are  transformed,  respect¬ 
ively,  to 

'  KVfl-F  + 1)  =  k*  .  m  =  m, ; 


I  These  relationships  yield  the  theoretical  relationships  between 

the  coefficient  of  friction  and  the  exponent  and  the  pressure  gradient, 
which  are  of  some  Interest.  Figure  1  shows  these  functions,  together 
with  the  experimental  points  pertaining  to  the  coefficient  of  friction 
j  .  [12]  In  the  following  coordinates;  c-/c-  ,  m/m^  and  a  (the  subscript 

I  ^ 

i  0  denotes  the  absence  of  a  gradient), 

i  b)  Let  f  >  0  and  a  =0  (injection  on  a  porous  plate).  Then  ac- 

I  '  ■ 

j  cording  to  (4),  (8),  (9),  and  (2) 

j  •  m  =  m^M  {.Y), 

I  e  _• 

I  -±  =  1(7)/,  ‘+'", 

}  .  ■  . 

where 

M(?)-2  I  inF|(l-rv. 

I  7'(F)  - 

•  The  . values  of  Y  are  found  from  the  following  equation; 

In  which  the  left-hand  part  Is  assumed  to  be  given. 


Pig.  1  Pig.  2 


The  calculations  that  have  been  carried  out  showed  that  despite 
the  absence  of  some  degree  of  freely  selected  constants,  the  experi¬ 
mental  values  of  m  In  the  most  varied  of  cases  are  In  satisfactory 
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agreement  with  the  theoretical  values,  provided  the  injection  velocity 
is  not  too  great  (at  high  injection  velocities  discontinuity  phenomena 
are  observed  and  the  corresponding  relationships  do  not  accurately 
describe  the  velocity  distribution);  an  analogous  comparison  of  frlc- 

f 

tlon-ooefflolent  values  is  difficult  In  view 
. of  the  substantial  scattering  of  the  experi-  ‘ 
mental  points.  As  an  example.  Pigs.  2  and  3. 
show  the  data  pertaining  to  Reference  [5];  the 
curve  in  Pig.  2  has  been  constructed  in  a  con¬ 
ventional  system  of  coordinates  c-/c-  and 

0 

2w/u  c~  (the  subscript  0  denotes  the  imper- 
0 

meable  surface)  and  corresponds  to  the  constant  Reynolds  ntimber  R^  = 

=  (UgX/v).  . 

c)  Let  f  >  0  and  a  <0  (injection  on  a  surface  streamlined  by 
an  accelerated  flow).*  In  this  case  the  sought  quantities  can  be  pre¬ 
sented  in  the  forro  of  m  =  m^MC?;  Z^,  Z^), 


where  M=,Z|(1— f  ^  K*/VI‘+" ,  Zi  =  K/*,  Z,  —  Ko*.  and  the  values  of 
Y,  Zj  and  Zg  are  determined  from  the  following  equations: 


P  dr  (  2  j 


and  Condition  (4)i  Solutions  taking  this  form  can  be  used,  for  exam¬ 
ple,  in  investigating  flow  in  the  vicinity  of  the  critical  point.  Plg- 
ure  4  shows  the  results  obtained-  in  the  calculation  of  the  frlctldn 
coefficient  for  various  local  values  of  the  pressure  gradient  in  the 

case  of  the  same  Reynolds  number  R  =  (u^S/v)  =.  10^  (the  subscript  0 

#  #■ 

denotes  the  zero  values  of  the,  parameters  f  and  a  ),  It  should  be‘ 
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pointed  out  that  the  values  of  the  coefficient  of  friction,  pertaining 
to  the  same  R  number  (for  example,  to  one  and  the  same  point  on  the  ' 
surface),  will  Increase  substantially  more  Intensely  with  an  Increase 
In  Sp/dx  than  shown  in  Pig.  4  as  a  result  of  the  reduction  In  the  R 

2.  The  outlined  concepts  can  be  used 
to  calculate  the  rather  general  case  of  a 
turbulent  boundary  layer  In  the  presence 
of  diffusion  and.  heat  transfer.  The.  ques¬ 
tion  Is  resolved  particularly  simply  If  we 
have  a  condition  of  similarity  between  the  , 
velocity,  temperature,  and  concentration 
^  fields,  1. e. ,  If  there  Is  a  "generalized" 

Crocco  Integral  (derived,  for  example.  In 
Reference  [8]).  In  the  general  case.  If  the  functions  governing  the 
changes  In  the  physical  properties  of  the  gas  are  known  (viscosity, 
thermal  conductivity,  etc.).  It  Is  always  possible  to  solve  the  system 
of  equations  which  describes  the  flow  near  the  wall  (the  equations  In 
Reference  [13]  evaluated  for  low  values  of  ^)j  for  this  It  Is  suffi¬ 
cient  to  pose  the  Koshl  [sic]  problem,  assuming  Inadequate  boundary 
conditions  at  the  surface  (for  example,  the  magnitudes  of  friction, 
heat  flow,  and  concentration).  The  solutions  found  with  the  aid  of  re¬ 
lationships  of  the  type  of  (4),  (8),  and  (9),  and  Assumption  (5)  ,  If 
It  actually  Is  valid,  make  It  possible  to  determine  the  wall  re-  ■  • 
glon  and  the  velocity  profile  In  the  turbulent  part.  With  respect  to 
the  determination  of  the  sought  magnitudes  on  the  surface  we  would 
point  out  that  the  equations  of  turbulent  flow  In  a  known  velocity- 
field  and  given  known  functions  governing  changes  In  "turbulent" 
Prandtl,  Schmldt>  etc.  numbers,  make  It  possible  to  associate  the  "ex¬ 


number. 


ternal"  boundary  conditions  with  the  conditions  prevailinjg  at  the 
boundary  of  the  viscous  region;  the  latter,  however,  in  view  of  the 
continuity  of  the  corresponding  magnitudes  (velocity,  temperature, 
concentration,  and  flows  of  momentum,  heat,  and  matter)  are  expressed 
In  terms  of  the  boundary  conditions  at  the  wall. 

As  an  example  let  us  consider  the  flow  of  a  binary  mixture  in  a 
boundary  layer  on  a  porous  surface  with  a  nonzero  negative  pressure 
gradient;  for  the  sake  of  simplicity  we  will  assume  that  there  is  no 
heat  transfer. *  .The  solution  for  the  viscous  region  in  this  case  will 


be  written  in  the  following  form: 


»;+/* 

j  '  /..Sc 


(10) 


where  y  =  i!^,. (p^  is  the  density  at  the  wall).  Sc  is 
the  Schmidt  number,  —  ^  1,1-)- (n —  1)C*) (p,  =  is  the  parameter  of 


the  pressure  gradient,  -4  =  (fi— i) 


I _ 

— — —  (C  and  C,  is  the  concentra- 


tlon  and  its  value  at  the  wall,  |x  =  M^/Mg  is  the  ratio  between  the 

molecular  weights  of  the  basic  and  introduced  gas).  In  the  derivation 

of  these  expressions  we  made  use  of  the  balance  equation  fpr  diffusion 

dC 


flows  PuWll  — C*)  =  —  p„D„ -^j  (^12  coefficient  of  diffusion); 

for  the  sake  of  simplicity  it  was  assumed  that  Sc  =  v/D^2  =  const  and 

as  a  result  as  well  as  in  view  of-  the  fact 
that  D-|^2  calculated  by  the  Enskog  [^]  method 

,1s  a  weak  function  of  the  mixture  proportion, 

* 

It  was  assumed  that  v  =  const. 

In  simplified  terms,  in  the  equations  of 
turbulent  motion  we  will  assume  d/dx  =  0 
.(for  example,  let  this  be  as  a  first  apprdxl- 


-  10 


I 


matlon);  then,  using  the  condition  at  the  boundary  of  the  regions,  we 
will  obtain 


—  dp  ^  ‘ 

+—  flf  — a 


mu^p 


+  D, 


(11) 


where 


“  *  =  (r?  — ^  +  - 


Se^  =  const  Is  the  turbulent  Schmidt  number. 

Integrating  the  second  of  the  equations  In  (ll),  we  will  find 
the  concentration,  distribution  which  as  a  result  of  the  condition  C 
=  0  in  the  case  of  y  =  6  yields  the  relationship 


n.-r 


PyW  J  ’ 


(12) 


where 


<D(/)-SCtJ  — 


d( 


dp 


dx  pyWU, 


For  a  given  distribution  of  w(x)  and  Ug(x)  Formulas  (4),  (8), 

(9),  (lO)j  and  (12)  make  it  possible  to  determine  all  of  the  unknown 
flow  parameters  contained  in  them.  The  results  of  the  calculations 
carried  out  In  the  assumption  that  dp/dx  =  0  and  Sc  =  Sc^  =  1  for  a 

mixture  of  air  and  helium  are  shown  in  Pig.  5.  All  of  the  values 

P  II 

,(c-/2  =  T  /p  u^)  pertain  to  the  same  Reynolds  number  R  =  10  .  Figure  5 

Illustrates  the  sharp  reduction  in  the  coefficient  of  friction  In  the 

case  of  the  Injection  of  a  light  foreign  gas;  however.  It  should  be 

pointed  out  that  the  final  results  of  the  relationship  between  the 
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physical  properties  of  the  mixture  and  the  concentration  (D2^2  “  const, 
M>  ,=  pScDj^2  =  P  const)  that  were  assumed  for  the  sake  of  simplicity  re¬ 
sult  in  a  somewhat  too  rapid  change  in  viscosityj  and  this  contradicts 
the  results  of  kinetic  theory.  Formally  this  indicates  that  with  the 
introduction  of  light  gases  (|x  >  l)  in  the  vicinity  of  y^^f  =  0  the 
derivative  d  u/dy  is  always  negative,  in  which  connection,  given 
small  flow  rates  for  the  ga,s  being  introduced,  there  is  observed  the 

unlikely . Increase  in  the  ratio  c„/c„  . 

I  '  ^  Q 

^  Received  10  July  I96I 

[Footnotes] 

Manu¬ 

script 

Page 

No. 

I  .  •  , 

1  In  Reference  [4]  it  has  been  suggested,  however,  that  the 

obtained  solutions  be  used  to  describe  the  flow  in  the  vi¬ 
cinity  of  the  critical  point. 

4  A  special  case  of  Relationship  (4)  was  used  in  the  work  of 

Van  Drelst  [4]. 

6  The  value  of  mQ  in  view  of  Assumption  (5)  need  not  be  a 

[  function  of  the  Reynolds  number,  and  for  this  reason  Formula 

(8)  does  not  take  into  account  the  Insignificant  evolution 
of  the  velocity  profile  in  the  case  of  constant  parameter 
values  (for  example,  on  the . Impermeable  plate). 

7  The  region  of  positive  pressure  gradients  is  not  considered, 
since  in  this  region,  generally  speaking,  the  assumption  of 

^  a  Single-parameter  velocity  profile  is  not  satisfied. 

10  Qualitative  evaluations  have  shown  that  in  accordance  with 

experimental  data  and  various  theoretical  works  the  effect 
of  heat  transfer  on  c^/c^^  and  St/Stg  (St  and  St^  are  the 

Stanton  numbers)  as  functions  of  the  corresponding  .dimen¬ 
sionless  mass  flow  rate  (see  Pig.  5)  are  generally  not  great. 
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THE  PROCESSES  OP  IONIZATION  AND  RECOMBINATION  IN  HQT  AIR 

V.  A.  Bronshten 
(Moscow) 

In  the  passage  of  extremely  powerful  shock  waves,  not  only  molec¬ 
ular  dissociation  can  take  place  in  the  air  behind  the  shock-wave 
front,  but  the  ionization  of  the  atoms.  An  example  of  such  a  shock 
wave  is  the  wave  formed  by  a  meteorite  flying  through  the  atmosphere 
(as  is  well  known,  the  velocity  of  meteorites  in  the  atmosphere  ranges 
between  11  and  72  km/sec). 

TO  calculate  the  processes  of  ionization  and  recombination  in  the 
hot  air  behind  the  front  of  a  wave,  the  air  in  a  plasma  state,  we  must 
know  the  coefficients  of  ionization  and  recombination  (r  is  the  ^ 
degree  of  ionization).  In  order  to  obtain  ^hese  coefficients  we  must 
evaluate  the  comparative  role  played  by  the  various  elementary  proc¬ 
esses,  and  in  contemporary  scientific  literature  this  is,  by  no  means, 
done  in  a  uniquely  defined  manner. 

1.  Ionization.  The  basic  mechanism  of  ionization  in  hot  air  is 
the  Impact  of  electrons.  The  general  expression  for  the  coefficient  of 
ionization  resulting  from  electron  Impact,  as  used  in  courses  in  astro¬ 
physics  [1,  2],  takes  the  following  form; 

*-■  0»  •• 

Here  v  is  the  velocity  of  the  colliding  electron,  Vq  Is  the  mini¬ 
mum  velocity  resulting  in  ionization,  k^^^  effective  cross 
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section  of  the  r-th  ionization  from  the  k-th  level,  at  velocity  v. 

In  this  case  it  is  generally  assumed  that  the 
ionization  takes  place  primarily  from  the  base 
level,  since  the  population  of  excited  levels 
In  the  case  of  the  Boltzmann  distribution  is 
negligibly  small.  Let  us  try  to  calculate  the 
from  the  base  level.  Generally  [1,  2] 

is 


-r,  1 


Pig.  1 

coefficient  of  ionization  Z 

the  concept  of  the  mean  effective  ionization  cross  section  "a. 
used  to  calculate  the  integral  in  (l),  placing  this  mean  value  outside 
of  the  radical  sign.  Then 


r,  V 


(2) 


2 

Here  it  has  been  taken  into  consideration  that  =  (m^VQ/2). 
However,  0  is  a  function  of  the  temperature  T  and  the  ionization 

X  V  6 

potential  1^^,  i.  e.  ,  a  function  of  v^.  It  is  therefore  more  convenient., 
as  was  ■  done  by  S.'B.  Pinkel'ner  [3]^  to  make  use  of  the  circumstance  that 

in  the  energy  Integral  of  Interest  to  us  the  quantity  .0^  ^  increases 

2 

almost  linearly  in  the  v  function  (Pig.  l). 

In  this  case 


a. 


(3) 


where  Oq  Is  the  effective  ionization  cross  section  at  an  energy  of 
2Ip,  i. e. ,  at  a  velocity  of  Vq.^2.  Substituting  (3)  into  (l)  and  cal¬ 
culating  the  Integral,  we  will  find  the  following  expression  for  the 
Ionization  coefficient.  * 


('*) 


Prom  a  comparison  of  (2)  and  (4)  it  follows  directly  that 
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1+ 


2kT, 


1+ 


n,  l,  +  2kT, 
If  + 

*r. 


(5) 


Substituting  the  numerical  values  of  the  constants  into  Formula 
(4)>  we  will  bring  it  to  a  form  convenient  for  calculations; 

If 

(4a) 


The  values  of  the  coefficients 
of  ionization  from  the  base  level, 
calculated  for  air  according  to  For¬ 
mula  (4a)  aj?e  presented  in  Fig.  2.  * 
However,  the  ionization  from 
the  base  level  by  no  means  determines 
the  general  course  of  the  ionization. 

As  has  been  shown  by  the  Works  of  G. S. 

Fig*  .2  Ivanov- Kholodnyy,  G. M.  Nlkol'skly,  . 

and  R. A.  Gulyayev  [4,  5],  as  well  as  in  the  works  of  L.M.  Biberman  and 
his  coworkers  [6,  7],  the  chief  contribution  to  ionization  is  made  by 
the  upper  levels.  Indeed,  despite  the  fact  that  the.  population  of  these 
upper  levels  is  small  in  comparison  with  the  base  level,  the  coeffi¬ 
cient  of  ionization  from  the  upper  levels  Increases  sharply  as  a  result 
of  the  reduction  in  the  binding  energy  of  the  levels. 

In  particular,  the  probability  of  ionizing  hydrogen  from  the  n 
level  is  equal  [6]  to 

where  u^  »  (ip/lcTn^) ,  and  -El(-Uj^)  is  the  integral  exponential  func¬ 
tion.  The  shape  of  the  product  of.  .  with  the  number  of  the 

n  level  for  various  values  of  u^  =  (l^/kT)  is  presented  |.n  Fig.  3. 
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»  tim  f 


To  determine  the  over-all  coefficient 
of  hydrogen  Ionization  (from  all  levels) 
is  summed  in  accordance  with  the  Boltzman 
distribution: 


i 


(/(n 


(7) 


where  is  the  statistical  weight,  U(t)  Is 
the  sum  of  the  given  atom  or  ion  according 


TABLE  1 


1  Bna  iMaSMU 

T 

»• 

Bm  njiasMM 
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r 

^PaaoBbie  ryMaHHoctH 
3MoKoct|wpa 

40MiHe«iHaii  KopoHa 

3  XpoMO€<|«pa 

10* 

10« 

l(f 

10»» 

1 

10* 

10» 

10* 

5-10» 

870 

420 

190 

58 

nporytfepaHufci  '"6 

ropRNHA  BosAyx  7 

To  )Ke  8 

III 

7,5- 10» 
10* 
to* 

40 

8—10 

1—2 

l)  Type  of  plasma;  2)  gas  nebulae;  3)  ionosphere;  4) 
solar  corona;  5)  chromosphere;  6)  protuberances;  7) 
hot  air;  8)  the  same. 


to  state.  Sum  (7)  exhibits  a  divergence  if  we  do  not  take  Into  consid¬ 
eration  the  deviation  of  the  distribution  over  the  levels  from  the 
Boltzmann  distribution,  as  a  result  of  the  mutual  excitation  of  the 
Ions  (the  effects  of  preionlzatlon  and  supercharging).  This  phenomenon 
was  examined  in  the  classical  Tomson  [sic]  approximation  by  G. S. 
Ivanov-Kholodnory  et  al.  [4]  and  In  the  Bethe-Born  approximation  by  L. 
M.  Biberman,  Yu. N.  Toropnlk  and  K.  N.  Ul'yanov  [6].  The  qualitative 
results  of  these  two  works  coincide  and  can  be  reduced  to  the  fact 
that  levels  above  sum  Uq  in  a  real  gas  cannot  be  achieved.  However,  an 
evaluation  of  n^  for  hydrogen  in  Reference  [4]  is  somewhat  underesti¬ 
mated  (i.e. ,  the  effect  of  supercharging  Is  exaggerated),  as  Is  shown 
by  a  comparison  against  experimental  data  [6], 

For  purposes,  of  evaluating  nQ  Reference  [4]  proposes  the  follow- 
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ing  formula  (close  to  the  analogous  Unzol'd  [sic]  formula) 

log  n,  =21,65— 6  log  (Hq  +  1), 


where  Is  the  electron  concentration.  For  the  various  types  of  plas¬ 
ma  In  [4]  the  following  values  of  Uq  have  been  found,  and  to  these  we 
have  added  the  values  for  hot  air  (Table  l). 

As  can  be  seen,  a  comparatively  small  quantity  of  levels  can  ac¬ 
tually  be  achieved  In  hot  air.  On  the  one  hand  this  circumstance  fa¬ 
cilitates  the  summing  and  makes  it  possible  to  do  entirely  without  the 
Integration  of  Formula  (7)  with  respect  to  the  upper  levels,  as  Is 
generally  done.  However,  the  small  number  of  levels  Imposes  greater 
requirements  upon  the  correct  evaluation  of  n^. 

L. M.  Blberman  and  his  coworkers  [6]  obtained  the  following  ap¬ 
proximate  expression  for  (in. the  Interval  from  0.01  <  u^  <  5,  where 
2 

Uo  =  I^/kTnQ)  as  a  result  of  Integrating  Formula  (6)  with  respect  to  ■ 
the  Boltzmann  distribution , (O  <  n  <  nQ) : 


(9) 


which  cam  be  transformed  to 


(10) 


In  the  astrophyslcal  examples  collected  In  Table  1  am  error  of 
2-3  units  In  the  evaluation  of  the  n^  number  has  little  effect  on  the 
result,  since  n^^  Is  great.  However,  In  our  case  when  n^  Is  small  an 
error  of  even  a  single  unit  results  In  a  change  In  by  an  entire 
order  of  magnitude,  since  nQ  enters  the  formula  raised  to  the  fifth 
power.  At  the  same  time.  In  a  real  gas  Hq  levels  will  never  be  achieved 
for  all  gas  atoms  simultaneously.  There  will  be  a  certain  statistical 
distribution  of  atoms  over  nQ  and,  generally  speaking,  the  effective 
value  of  nQ  will  be  expressed  by  a  fraction.  Below  we  discuss  the  me- 
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thod  of  evaluating  it. 


2.  Recombination.  Now  let  us  examine  the  phenomenon  of  recombina¬ 
tion.  The  latter  may  take  place  in  one  of  two  ways;  with  the  emission 
of  a  radiation  quantum  and  as  a  result  of  triple  collisions,  in  which 
case  the  third  particle  which  transmits  the  excess  energy  will  almost 
always  be  an  electron. 


t 


I 


For  the  coefficient  of  recombination  with  radiation  in  the  hydros 
gen-similar  approximation  we  can  use  the  following  well-known  formu¬ 
la  [8,  91:* 


/  2*r,  X  L 

’  3!L(3)* 


(11) 


where  z  is  the  ion  charge  (in  proton -charge  units),  is  the  capture 
constant: 


16 

3/3  mjc* 


=.2^11.10-M  CM*; 


9(p)  is  the  function  presented  in  tabular  form  by  Spitser  [9J  and 
equal  to; 


(12) 


The  quantity’ <p(p)  for  the  required  interval  T^  is  of  the  order  of 
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p 

unity,  as  can  be  seen  from  Pig.  4.  Evidently,  p/n  =  u^. 

The  comparative  role  of  the  various  levels  Is  here  determined  by 
the  temperature.  With  an  Increase  In  temperature  the  relative  quantity 
of  recombinations  at  the  base  level  will  increase  [4];.,  Therefore  the 
above-mentioned  restriction  of  the  number  of  levels  has  little  effect 
on  recombination  with  radiation  (see  the  curves'  In  Pig.  5  for  the  co¬ 
efficients  of  recombination  with  radiation  for  air).* 


Por  the  third  recombination  (or  recombination  with  triple  colli¬ 
sions)  In  Reference  [6]  there  Is  derived  the  formula  which  determines 
the  probability  of  recombination  at  the  n  level: 

and  after  Integration'  with  respect  to  the  levels  we  obtain 

or,  after  substitution  of  (lO)  Into  (l4): 

\  kT  f 

However,  we  must.  In  this  case,  take  Into  consideration  that  a 
portion  of  th©  upper  levels  Is  not  attained  as  a  result  of  oncoming 
Ionization,  and  therefore  the  effective  coefficient  of  triple  recom¬ 
bination  takes  the  following  form: 


c(») 

c* «» y 

'  f  ‘  +  ^r^V4 


(16) 


where  is  the  probability  of  spontaneous  transition  from  the  n 


level  (A  =  1.6*10^°  h”^'  . 

The  authors  'of  Reference  [6]  suggest  the  summing  to  some  level 
n^  <  beginning  from  which  the  terms  in  Sum  (l6)  sharply  diminish, 
ajid  the  units  in  the  denominator  can  be  neglected.  The  remainder  Is 
determined  by  means  of  integration  from  n^  to  nQ. 

For  the  purposes  of  our  problem,  however,  this  method  is  not 
suitable  because  of  the  considerations  described  above.  Therefore  we 
will  employ  the  concept  of  the  effective  boundary  level  with  the  quan¬ 
tum  number  n^^  which  we  will  define  as  the  quantity,  which  after  being 
substituted  into  Fomiula  (15)  in  the  place  of  Uq^wIII  yield  the  exact 
value  of  C  calculated  according  to  Pornrula  (l6)  by  direct  summation. 

As  an  example  we  will  sum  the  values  of  in  the  case  of  quadru¬ 
ple  air  ionization  (l^  =  77  ev)  at  T^  =  10^°K, • assuming  for  the  time 
being  that  all  of  the  levels  are  "hydrogen-similar"  (this  is  necessary 
for  purposes  of  Illustrating  the  method). 

As  we  can  see  from  Table  2  the  effect  of  oncoming  ionization  re¬ 
sulting  in  the  disruption  of  the  upper  levels,  is  particularly  pro- 
nounced  in  the  case  of  great  n  .  The  maximum  value  of  C  ,  with  an  in- 

crease  in  n^,  transfers  from  the  third  level  to  the  second,  and  then 

20 

to  the  base  level.  The  value  of  n^^  for  n^  >10  is  close  to  unity 
or,  in  other  words,  the  "lifespan"  of  the  upper  levels  is  rather  small. 

Figure  6  shows  the  effective  coefficient  of  triple  recombination. 
C^^^ng  (solid  line)  and  n^^  (dashed  lines)  to  be  functions  of  the  elec¬ 
tron  concentration  n  (l  r-  hydrogen,  T„  =  17,500°,  I_  =  13. 56  ev;  2  -  . 

^  I* 

air,  T  =  100,000°,  I  =  77  ev,  quadruple  Ionization).  The  ratio  I„/kT- 
is  equal  for  both  of  the  examples.  We  obtain  the  following  pattern 
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(pig.  6).  In  the  case  of  rather  small  concentrations  the  quantity 

(3) 

C'  'n^  Increases  slowly  —  there  is  a  transition  of  the  recombination 
maximum  from  the  upper  levels  downward,  to  the  base  level.  Then  there 
is  a  rapid  rise  in  C'  'n^,  since  despite  the  weakening  of  the  role 
played  by  the  upper  levels  the  oncoming  ionization  from  the  base  level 
still  does  not  play  any  peirticular  role  and  C^^^  remains  almost  con¬ 
stant  as  n  increases.  Finally,  with  a  further  Increase  in  n  a  sata- 
ration  sector  sets  in,  and  this  can  be  attributed  to  the  oncoming 
ionization.  In  this  same  Fig.  6  is  also  shown  the  shape  of  n^^  which 
approaches  unity.  With  a  drop  in  T,  n^^  Increases  approximately  as 
1-1/7. 


TABIE  2 


♦  ' 

rt,  »I0'« 

1 

4,2*10-» 

\8.4.10‘»»  ! 

8,4-10-« 

8,4.10-" 

8,3.10-" 

2 

2*4. 10-*  . 

2,3’i0“>« 

2,2’10"»* 

5,1.10-»> 

6,5.10-” 

3 

2,3-10-»’ 

1.4.10-” 

4,7.10-*» 

7,0.10-« 

7,0.10-” 

4 

7,0- 10-* 

4,9.10-«* 

2,1.10-»* 

2.2.10-” 

2,2.10-” 

5 

1,4  10-’ 

1,2- 10-“ 

1,0- 10-” 

1,0.10-” 

6 

2,2.10-’ 

2,6- 10-" 

5,8.10-« 

5,8*10’« 

5,8.10-” 

l,3  10-’» 

3,23 

7.2.10-« 

1,86 

9,1.10-" 

1.24 

The  comparative  role  of  recombination  with  radiation  and  triple 
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recombination  depends  on  the  value  of  n^  and  the  magnitude  of  ioniza¬ 
tion.  For  example,  at  '=  10^^  for  the  Ions  NV  and  OV  the  triple  re- 
combination  is  1-2^  of  the  recombination  with  radiation  and  It  can,  be 
neglected.  At  n  =  lO"*"  It  amounts  to  and  at  n  =  10  it 

amounts  to  120$^  of  the  recombination  with  radiation.  With  transition 
to  the  lower  ions  the  role  of  triple  recombination  Increases. 

3.  Transition  to  complex  atoms.  Up  to  this  point  we  have  regarded 
all  levels  in  the  atoms  and  ions  to  be  "hydrogen -similar"  although  in 
actual  fact  this  is  by  no  means  always  the  case.  In  order  to  make  the. 

calculation  for  complex  atoms  and  ions  more  exact  we  can  use  the  con- 

* 

cept  of  the  effective  quantum  number  n]^,  determined  from  the  following 
condition. [8,  10]: 


where  is  the  binding  energy  of  the  level  defined  by  the  quantum 
numbers  n  and  1;  I„  =  13. 56  ev  represents  the  potential  of  ionization 

“  In  ,  . 

for  the  hydrogen  atom;  z  =  r  +  1  is  the  charge  of  the  atomic  remainder 
(z  =  1  for  neutral  atoms,  z  =  2  for  single  ions,  etc.).  The  difference 
n  —  n  =11,  referred  to  as  the  quantum  defect,  vanishes  with  an  in¬ 
crease  in  levels,  converting  the  upper  levels  into  "hydrogen -similar" 
levels. 

In  the  works  of  L.M.  Blberman  and  his  coworkers  [6,  7>  12] 

there  are  developed  the  formulas  by  means  of  which  we  can  take  Into 
consideration  the  quantum  defect  for  the  actual  levels  of  complex 
atoms,  and  these  formulas  also  make  it  possible  to  find  for  these 
atoms  the  coefficients  of  recombination  and  Ionization.  The  work  by 
Burgess  and  Seaton  [10]  presents  the  general  formula  for  the  photo- 
lonlzatlon  section  of  complex  atoms  and  hvimerous  examples  of  how  this 
formula  can  be  used. 
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Pig.  8 

The  formula  for  the  coefficient  of  recombination  with  radiation 
for  complex  atoms  takes  the  following  form; 


C,  =1.04  1(f"r  ’  (  -t(yr  ^).  e-"(e"'—iUv, 
J  M 


(18) 


viiere  the  frequency  v  is  determined  by  the  energy  of  the  recombining 

electron  (hv  =  mv  /2),  u'  =  hv'/kT,  v'  =  v,  if  v  <  v_  and  v '  =  v  ,  if 

s  s 

V  >  V  ,  where  v  is  the  frequency  corresponding  to  the  lower  excited 
g  s 

state  n  from  which  the  integration  begins  (regarding  the  selection 

s 

of  n  and  v  see  [11]). 

o  o 

The  difference  from  "hydrogen -similarity"  is  expressed  by  the 
function  4(v,  T) ,  which  for  various  atoms  and  ions  exhibits  a  varied 
shape  [11].  For  atoms  of-  oxygen  and  nitrogen  ^(v)  <  1  is  not  a  func¬ 
tion  of  temperature  and  diminishes  almost  linearly  With  v(for  v  < 


.-15 


<  lO"'^-^  sec“^).  But  for  various  ions  of  oxygen  the  shape  of  ^(v)  va¬ 
ries  substantially,  and  in  this  case  this  quantity  may  be  more  or  less 
than  unity  [12].  However,  with  an  increase  in  the  multiplicity  of  the 
ions  the  shape  of  i(v)  is  smoothed,  approaching  the  straight  line  4  = 
=1,  and  this  is  explained  by  the  fact  that  the  levels  of  ions  of  high 
multiplicity  approach  ever  closer  to  the  "hydrogen-similar"  levels. 
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At  the  high  temperatures  taking  place  in  the  air  in  the  case  of 
shock  waves  in  front  of  a  flying  meteorite  [T  =  (5-20)* 10^"K],  the 
first  ionizations  take  place  quickly  and  do  not  make  a  great  contribu¬ 
tion  to  the  specific  heat,  whereas  the  recombination  of  ions  of  small 
multiplicity  is  negligibly  small.  Therefore,  in  this  problem  the  uti¬ 
lization  of  the  exact  Formula  (l8)  is  not  Justified  by  need,  and  we 
can  assume  everywhere  that  4=1.  The  possible  errors  for  the  OI-OIII 
ions  partially  compensate  each  other,  since  the  difference  4—1  for 
these  ions  is  opposite  in  sign. 

For  our  calculations  we  will,  therefore,  employ  Fonnulas  (ll)  and 

(12). 

With  respect  to  ionization  and  triple  recombination,  the  applica¬ 
tion  of  Formulas  (lO)  and  (l6)  to  the  complex  atoms  need  also  not  re¬ 
sult  in  large  errors  [6].  However,  the  factor  Tj,  must  be  Introduced 
into  Formula  (lO),  since  this  factor  makes  it  possible  to  take  into 
consideration  the  multiplicity  of  the  atoms ; 


where  is  the  s\im  over  the  states. 

Figures  7  and  8  show  the  values,  respectively,  of  the  coeffi¬ 
cients  of  ionization  from  all  levels  and  of  the  coefficients  of  triple 
recombination  for  air,  calculated  according  to  the  Indicated  formulas 
for  the  electron -temperature  Interval  T^  =  (3-30)  19^‘’K. 

The  author  wishes  to  express  his  gratitude  to  K. P.  Stanyukovlch 
for  his  continuous  attention  to  the  work,  and  to  L.M.  Blberman  and  G. E. 
Norman  for  their  useful  discussions  and  the  fact  that  they  made  it 
possible  for  the  author  to  familiarize  himself  with  materials  con¬ 
tained  in  the  manuscripts  of  articles  [6,  7);  In  addition,  the  author 
wishes  to  express  his  thanks  to  S.B.  Plkel'ner  for  his  valued  advice. ' 
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[Footnotes] 
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16  In  Reference,  [3]j  In  the  place  of  the  factor  (1+2  (kT^/l^)) 

there  Is  simply  and  this  Is,  apparently,  an  error. 

.Therefore  the  coefficients  are  underestimated  in  this  work 

approximately  0. 5  orders  of  magnitude,  which,  by  the  way,  has 
little  effect  on  the  results  of  this  work  because  of  the  pre¬ 
dominant  role  played  by  the  electron  excitation  In  nebulae. 

17  In  all  calculations  air  was  assumed  to  be  a  monocomponent 
gas,  and  the  potentials  and  coefficients  of  Ionization  for 
nitrogen  were  averaged  In  accordance  with  the  relative  con¬ 
tents  of  these  gases,  l.e. ,  0. 78  and  0.22,  respectively. 

20  In  the  work  by  Seaton  [8]  this  formula  has  a  slightly  dif¬ 
ferent  form,  but  it  Is  hot  difficult  to  modify  it  to  the 
form  of  (11).  However,  this  formula  Is  suitable,  strictly 
speaking,  only  for  the  Isoelectron  hydrogen  series. 

21  The  question  of  recombination  with  radiation  for  oxygen  and 
‘nitrogen  Is  also  discussed  In  Reference  [13]. 
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AERODYNAMIC  CHARACTERISTICS  OP  DELTA  WING  WITH  DEFLECTED  AILERON, 
FLAP,  AND  OTHER  SECONDARY  CONTROL  SURFACES  IN  SUPERSONIC  FLOW 

Yu. I.  Krasil'nikov 
(Moscow) 

In  this  work  we  obtain  the  aerodynamic  characteristics  of  a  delta 
wing  with  a  deflected  aileron,  flap,  and  other  secondary  control  sur¬ 
face,  when  streamlined'  by  a  supersonic  flow  of  gas. 


The  theoretical  determination  of  the  aerodynamic  coefficients  Is 
based  on  the  linear  theory  of  a  wing  of  finite  span  In  a  supersonic 
stream,  as  developed  In  Reference  [1]. 

1.  Statement  of  the  problem.  Let  us  consider  the  movement  of  a 
supersonic  gas  stream  past  a  thin  delta  wing  with  deflected  aileron, 
flap,  or  other  secondary  control  surface.  The  wing  plane  will  form  a 
small  angle  of  attack  a  with  the  direction  of  the  velocity  tJ  of  the 
approaching  stream.  The  aileron,  flap,  or  other  secondary  control  sur¬ 
face  Is  deflected  from  the  plane  of  the  wing  through  the  small  angle  6 
(plg.  l).  The  leading  edges  of  the  wing  are  supersonic. 
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Let  the  flow'  being  formed  satisfy  the  conditions  of  the  linear¬ 
ized  theory  [ij.  Given  the  indicated  streamlining  conditions  we  will 
determine  the  aerodynamic  characteristics  of  the  wing. 

Let  us  introduce  the  basic  rectangular  system  of  coordinates  xo'^^ 
oonneotedi  to  th®  wing,  and  the  system  of  characteristic  coordlnat®# 
xoyz  (Fig.  l).  The  connection  between  the  characteristic  coordinates 
and  the  basic  coordinates  is  accomplished  in  accordance  with  the  fol-  ' 


loinlng  formulas  known  from  [1] 


■~x  —  ky,  y^x-^ky,  t  =  kt. 


(1.1) 


where  k 


=  - 


1  =  cot  e;  M  >  1  is  the  Mach  niunber  of  the  approach-. 


Ing  stream  and  e  is  the  Mach  angle. 

The  velocity -disturbance  potential  <p  satisfies  the  wave  equation. 
In  accordance  with  the  linear  theory  for  the  wing,  the  boundary  condi¬ 
tions  for  the  potential  9  in  the  characteristic  coordinates ,  after 
linearization,  will  be  the  following; 


on  the  projection  of  the  basic  part  of  the  wing  onto  the  xoy 


plane 


i  —  t/d  tg  e; 


on  the  projection  of  the  aileron,  flap,  or  secondary  control  surface 


•  =  — y  («  + A)fge. 


(1.2.) 


Since  the  boundary  conditions  (1.2)  are  constant,  the  value  of 
the  potential  <p  on  the  projection  of  the  wing  (on  the  basis  of  the 
general  formula  for  the  determination  of  the  potential)  in  the  charac¬ 
teristic  coordinates  is  written  in  the  following  form; 


’  '  ^  sA.\) 


^  ru-8)(i>-n)  ’ 


(1.3) 


tg  =  tan;  ctg  =  cot 


where  S2(x,  y)  Is  the  integration  region  with  respect  to  the  projec¬ 
tion  of  the  basic  part  of  the  wing  falling  within  the  characteristic 
cone  whose  apex  is  situated  at  the  point  M(x,  y) ;  and  Sg(x,  y)  is  the 
integration  region  with  respect  to  the  aileron,  flap,  or  secondary 
control  surface  projection,  falling  within  the  same  cone.  . 


To  determine  the  forces  and  moments  acting  on  the  wing,  it  is 
necessary  to  know  the  distribution  of  the  pressure  differences  across 
the  wing.  In  accordance  with  the  Bernoulli  integral,  the  difference 
between  the  pressure  beneath  and  above  the  wing  is  associated  with  the 
velocity -disturbance  potential  <p  by  the  following  relationship: 

AP  =  -2pt;(<p,+.p^).  (1.4) 

where  p  is  the  density  of  the  gas  in  the  undisturbed  stream  and  <p^  and 
<p  are  the  derivatives  of  the  potential  cp  with  respect  to  the  charac- 

y 

terlstic  coordinates. 

Thus  in  order  to  resolve  the  posed  problem  it  is  necessary  to  de¬ 
termine  the  velocity -disturbance  potential  T^and  then  to  find  the  dis¬ 
tribution  of  the  pressure  difference  across  the  wing. 

2.  Determination  of  the  pressure  difference,  a)  Distribution  of 
the  pressure  difference  across  a  wing  with  an  aileron.  In  flow  past  a 
delta  wing  with  a  deflected  aileron  (Fig,  la)  two  oases  are  possible. 
In  the  first  case  (Pig.  2a)  the  aileron  falls  completely  within  the 
head  Mach  cone,  whereas  in,  the  second  case  (Pig.  2b)  the  aileron  pro- 
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trudes  from  the  head  cone.  In  both  cases  the  disturbance  (perturba¬ 
tion)  lines  emanating  from  points  0,  0^,  and  0^  divide  the  entire  flow 
region  on  the  wing  into  regions  I,  II,  III,  IV,  V,  and  VI,  and  flow  in 
these  regions  differs  in  nature.  The  equations  of  the  wing  and  aileron 
generatrices  are  assumed  to  be  in  characteristic  coordinates  when  the 
potential  is  sought  in  these  regions. 

The  flows  in  regions  I  and  II  are  known.  The  distribution  of  the 
pressure  difference  in  these  regions  is  determined  on  the  basis  of  the 
following  formulas; 

AP,  =  4(7algef  (t.  e), 

=  Aqa.  (g  tF  (r,  e)  (t,  e,  <p). 

Where  -we  have  denoted 

<p  is  the  polar  angle  of  the  point  in  region  II  counted  off  from  the 
ox  axis,  and  here  <p  >  0  in  the  region  Oij^OE  and  <p  <  0  in  the  region 
O^OP;  c[  is  the  .velocity  head. 

Using  Formulas  (1.3)  and  (1.4),  we  obtain  the  following  expres¬ 
sion  for  the  pressure  difference  in  region  III: 

A^lll  =  ef  (t,  e)  <P,  (t,  e,  <p),  (2.5) 

where  9  is  the  polar  angle  of  the  point  in  region  III  counted  off' from 
the  ox  axis. 

In  seeking  the  potential  q)j.y.(x,  y)  in  region  IV,  we  do  not  take 
into  consideration  the  effect  of  the  slit,  but  we  do  take  into  con¬ 
sideration  the  leakage  of  the  gas  from  one  half  of  the  region  to  the 
other.  Integration  in  the  region  IV  is  more  conveniently  carried  out 
in  characteristic  coordinates  whose  origin  is  situated  at  the  point 


(2.1) 

(2.2) 

(2.  3) 
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For  the  pressure  difference  In  region  IV  we  obtain  the  following 
formula; 

APiy  =  <9®  Ig  e  -L  arc  tg  l/" -li-  + 
n  **  r  yt 


where 


(2.6) 


tg  Y  —  tg  e  I 

tgT  +  tge  m 


(2.  7) 


and  are  the  characteristic  coordinates  of  the  point  (the  origin 
Is  at  point  0^)  In  region  IV;  B,  b,  and  h  are  the  geometric  parameters 
shown  In  Pig.  1. 


Pig.  3 


The  connection  between  the 
characteristic  coordinates  (x,  y) 
whose  origin  Is  situated  at  point  0 
and  the  characteristic  coordinates 
(Xi,  y^^)  whose  origin  Is  situated  at 
point  Is  achieved  by  the  follow¬ 
ing  formulas ; 

jr  =  x,+;(fl  — 6)— ft®.  ».5=gi  +  (B-®)  +  *®-  (2.8) 


On  the  basis  of  (2.8)  Formula  (2.6)  can  be  rewritten  In  the  fol¬ 
lowing  form: 

AP,y  «=  <9®  »>■«  tg  + 

-f  (t. 

Rearranging  and  Introducing  the  denotation 


1- • 


(2.9) 


(2.  10) 


for  the  pressure  difference  In  region  IV  we  will  obtain  the  final  ex¬ 
pression 

AP,y  =  4gA  tg  etPi  (e,  (p,)  +  4<7»  tg  e®,  (t,  e,  y),  (2.  11 ) 
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where  9^^  Is  the  polar  angle  of  the  point  In  region  IV,  counted  off 
from  O^D  (the  pole  is  situated  at  point  O-j^) ,  and  9^  >  0  In  the  region 
O^DG;  9^  <  0  in  the  region  O^DKj  9  is  the  polar  angle  of  the  same 
point  In  region  IV,  counted  off  from  the  ox  axis  (the  pole  is  situ¬ 
ated  at  point  0).  For  the  remaining  regions  the  pressure  difference  Is 
■determined  by  the  following  formulas: 

^  APy  4- <p,).  (2.12) 

APy,  U*f  Ct.  e)  +  Mli!e  (2.13) 

taking  into  consideration  (2.3)  and  (2.4). 

b)  Distribution  of  pressure  difference  across  a  wing  with  flap 
and  secondary  control  surface.  In  the  case  of  flow  past  a  delta  wing 
with  a  flap  or  secondary  control  surface  the  entire  flow  region  Is 
divided  into  eight  characteristic  regions  I,  ll.  III,  IV,  V,  VI,  VEI, 
and  VIII  by  perturbation  (disturbance)  lines  emanating  from  points  0, 
Oj.,  and  Og  (pig.  3)  in  the  general  case. 

In  regions  I  and  11  of  a  wing  with  a  flap,  the  pressure  differ¬ 
ence  Is  determined  according  to  Formulas  (2. l)  and  (2.2),  taking  Into 
consideration  (2.3)  and  (2.4).  For  a  wing  with  a  secondary  control 
surface  we  will  have 

AP,  =  4?  (a -r- ahb  Bp  (t.  «)■ 

APii  =  ‘*9  («  4- 'R  *P  (T.  «)'< 

,  e,  (p). 

The  pressure  difference  In  region  VTI  of  a  wing  with  a  flap  is  de 
termlned  according  to  Formula  (2. l4)  and  according  to  Formula  (2.1)  in 
the  case  of  a  wing  with  a  secondary  control  surface.  In  region  III  of 
a  wing  with  a  flap  the  pressure  difference  is  determined  on  the  basis 
of  Formula  (2.5). 

For  a  wing  with  a  secondary  control  surface,  we  will  obtain 
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(2.14) 

(2.15) 


(2.  16) 


.  .  .  ' 

AP,„  =  4?  (o  +  t>)  tg  er(t.  e)  (I>,  (t.  *.  f)  —  *8  «•. 

In  region  VI  of  a  wing  with  a  flap  the  pressure  difference  Is  de¬ 
termined  according  to  Formula  (2.13)-  For  a  wing  with  a  secondary  con¬ 
trol  surface  we  will  have 

'  APvi» 4? («  -(-«)  fgef(T,e)-Vtge.,  (2.1?) 


For  the  pressure  difference  In  the  region  V  of  a  wing  with  a  flap 
we  obtain  the  expression 

APv=V<R*f(T.  *)  +  V »gr.(«n,(e.  <(),)  + P(T.e)Jt-<I»,(T.e.<pt)|)  (2.18) 


where 


I'j  (T.  e.  <Pi) 


==  JL  *fc  tg  l/ <*gT  +  »g  e)(tge  — tg«Pi) 
«  y  (IgT  — «ge).(tge  + tg<Pi) 


(2.19) 


Is  the  polar  angle  of  the  point  In  region  V  (the  pole  Is  situated 
at  point  0^).  For  a  wing  with  a  secondary  control  surface  we  will  have 
APv'=4g(«H-A)tg«P(T.e)  — 4gAfge{<l*,(f,<p,)T-  P(t.  e)(l—<t>j(T.  K. 

For  the  pressure  differences  in  regions  IV  and  VIII  of'  a  wing 
with  a  flap  we  obtain  the  following  formulas: 

/\P,v  =  iqa  Ig  eP  (t.  e)  <D,  (f,  «,  (p)  + 

-f4g61ge  {<Pi(e.  ip,)-t-P(T,  e)(t  —  e,  (p,)|).  (2.2l) 

APyiii  —  tg  ef  (t,  e)  ®,  (y,  e.  qi)3-t-’4vA  tg  e  t®.  (e,  (p,)  P  (t,  e)  {1  —  ®,  (t,  e.  qi,))  -|- 

(2.22) 

+ 1®,  (e.  <p,)  — 1( -I- P  (T,  e)  (1  —  ®,(r.  e.  <p,)I). 


where  <p  Is  the  polar  angle  of  the  point  In  region  VIII,  counted  off 
from  the  ox  axis  (the  pole  Is  situated  at  point  o),  9^^  is  the  polar 

angle  of  the  same  point  In  region  VIII,  counted  off  from  the  line  O^D 

.(the  pole  Is  situated  at  point  0^),  tpg  polar  angle  of  the  same 

point  In  region  VIII,  counted  off  from  the  line  O2C  (the  pole  Is  situ¬ 

ated  at  point  Og),.  and  92  >  0  In  the  region  OgCM  and  9g  <  0  In  the  re¬ 
gion  OgCN. 

In  order  to  derive  the  formulas  for  the  pressure  difference  in 
regions  IV  and  VIII  of  a  wing  with  a  secondary  control  surface,  it  Is  , 
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sufficient  to  replace  cx  by  (a  +  6)  and  6  by  —6  in  Formulas  (2. 21)  and 
(2.22),  which  follows  from  an  analysis  of  the  formulas  derived  above. 

3.  Determination  of  forces  and  moments  acting  on  the  wing.  The 
forces  and  moments  acting  on  a  wing  are  determined  by  integration  with 
respect  to  all  regions  on  the  wing,  and  this  is  most  conveniently  car¬ 
ried  out  in  polar  coordinates. 

a)  Characteristics  of  a  wing  with  an  aileron.  Since  in  all  of  the 
forrouTas  for  the  determination  of  the  pressure  difference  in  various 
regions  of  a  wing  with  an  aileron  there  are  terms  that  are  functions 
only  of  a,  we  will  separately  calculate  the  force  acting  on  the  wing 
and  determined  by  these  terms.  In  either  case  (Pig.  2a  and  b)  we  will 
.  obtain 

4<7a  tg  eB»  tg  T-  (3*l) 

Using  Fonnulas  (2. 5)  sind  (2. 13)  we  will  find  the  force  acting  on 
the  region  0^0202^K  and  dependent  on  6.  Carrying  out  the  calOulation,. 
we  will  obtain 

z;  =  /w«l«.e(6fc - (3.2) 

The  exact  same  force  acts  on  region  OgO^OiiN.  On  the  basis  of 
(2.11)  and  (2.12).,  with  (2.10)  taken  into  consideration,  we  will  de- 
termine  the  force  acting  on  region  K0-|^G  and  dependent  on  6.,  Carrying 
out  the  Integration,  we  will  find 

z'  =  4</atge-^tge.  (3*3) 

The  force  acting  on  the  region  GOj^OgM  and  dependent  on  6  is  de¬ 
termined  from  the  following  formula: 

Z,=  =  4gfttge24A.  ,  (3*  4) 

For  the  total  force  acting  on  the  wing  we  will  obtain  the, follow¬ 
ing  expression: 
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Z  tg  eB*  tg  7  -r  ®  26/i. 


(3.5) 


The  coefficient  for  this  force  has  the  following  form: 

c,yAi«-i  /.  ■■ 

a  r  /’ 


(3.6) 


'  2 

where  =  2bh  is  the  total  aileron  area,  =  B  tan  y  Is  the  total 
area  of  the  wing, 

Let  us  determine  the  moment  of  forces  with  respect  to  the  oy 
axis.  The  moment  of  forces  (3*1)  is  determined  in  accordance  with  the 
following  formula: 


M.  ^Zi 


(3.7) 


For  the  moment  of  forces  (3.^)  we  obtain  the  following  expres¬ 


sion: 


(3.8) 


The  total  moment  of  all  forces  with  respect  to  the  oy  axis  takes 
the  following  form: 

Al-  =  4i/atge  T  + 417a tge2flt/i  1^1  — (3.9) 


For  the  coefficient  of  this  moment  we  obtain  the  following  for¬ 


mula: 


(3.10) 


Dividing  (3. 10)  by  (3.6)-,  we  will  obtain  the  dimensionless  coord¬ 
inate  for  the  center  of  pressure  of  the  central  chord.  Carrying  out 
this  division,  we  will  obtain. 


''HA 


'  3 


(3.11) 


where 
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(3.12) 


Ajr„ 


±  (1.  JJ\  h. 

a  \  3  ~  2b  )  ' 


Let  us  determine  the  moment  of  forces  acting  on  one  half  of  the, 
wing,  about  the  ox  axis.  The  total  moment  of  forces  acting  on  this 
half  of  the  wing,  about  the  ox  axis,  takes  the  following  form 

tgT 


arc  cos 


|/tg*T  — tg*e 


tgT  J^. 


(3- 13)' 


For  the  coefficient  of  this  moment  we.  will  obtain  the  following 
expression: 

tg.Y  ‘gel, 

- arC  COS - 

(3.14) 


«_yA»nrr_ 

3  n  [  tg  V 


4a 


^  y  - — -  div  vw 

V‘g’T  — ‘g’e  ‘gT  J 

ft  1  /  1.  \  h  S. 

p 

Where  Sg  =?  bh  is  the  area  of  one  half  of  the  aileron,  =  (B  /2)  tan  -y 
Is  the  area  of  one  half  of  the  wing,  and  H  =  B  tan  y  Is  the  half -span 
of  the  wing';  X  =  b/h. 

Formula  (3.13)  can  be  used  for  the  determination  of  the  moment  of 
forces  of  a  rotating  wing  or  control  surface  about  the  side  chord  In 
the  case  of  flow  past  a  wall.  At  a  =  0  we  will  obtain  the  moment  pro¬ 
duced  by  the  deflection  of  the  aileron. 

The  dimensionless  coordinate  of  the  center  of  pressures  of  the 
forces  acting  on  the  half  of  the  wing,  about  the  half -span.  Is  found 
from  the  following  expression; 

(3.15) 


“  C,.' 


b)  Characteristic  of  a  wing  with  a  flap  or  secondary  control  sur¬ 
face.  Let  us  determine  the  forces  and  moments  acting  on  a  wing  with  a 
flap.  Using  the  formulas  for  the  pressure  difference,  obtained  for  a 


-  38  - 


wing  with  a  flap,  we  will,  through  integration,  find  the  force 


Z  =  Vtg.*fl»tg;T  +  4<fdlge(26A4-MtgtI.  (3-l6) 

For  the  coefficient  of  this  force  we  will  obtain  the  expression 

■  7,-  A  .  (3.17),  • 

p 

where  =  (2bh  +  b  tan  y)  is  the  area  of  the  flap. 

The  moment  of  force  (3.16)  about  the  oy  axis  takes  the  following 

form; 

•Af-  =  Ago.  tg  e  Y  BMg  T  4-  4?*  tg  e  {2bh  +  fci  tg  t)  ^  fl  —  ^  —  iqi  fg  e-^  b*h.  (  3 .  l8) 

For  the  coefficient  of  the  moment  (3. l8)  we  obtain  the  following 
formula; 

V'M  —  1  ^  I  ^  [  /  1  t  \  ^3  _  b  '■S,,n.'[ 

4a  -  3  a  3B  /  S,  6fl  S,  J-  (3.19) 


where 'S„’  =  2bh. 

z.  p 

Dividing  (3.19)  by  (3.17),  we  will  obtain 

_  2  - 
“  “3“^  ‘'•''n  -T'. 


where 


^x 


n.A 


A 

<  ±~\ 

S3 

b 

^3.n  1 

a 

3  ( 

B  J 

6B 

Sk  j 

(3.20) 


(3.21) 


The  total  moment  of  forces  acting  on  one  half  of  the  wing  about 
the  central  chord,  takes  the  following  form; 

tgr  '■ 


•'Al-=4gatge4-fl*tg*T— f  — +  ^ 


-arc  cos 


^7]- 


t4g43g.A{f-[l+.-^(-f]%,..j,+.i;^^ 


(3-  22) 


For  the  coefficient  of  this  moment  we  obtain  the  following  ex¬ 
pression: 


n-VM’-l 

X 

4ct 


^  J  2  r  tge. 
3  L  tgT 


tg  T  .  tge 

4-  — ==i====rarccos- - 
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,■+ 


where  we  denote 


+■- 


■i[i( 


Si  =  “'C  Y.  St  ==  hli.  .S, ^ 

The  dimensionless  coordinate  of  the  center  of  pressure  of’  the 
forces  actlhg  on  the  half  wing  about  the  half  •span  can  be  found  fro<n 
Formula  (3*15). 

In  order  to  determine  the  aerodynamic  coefficients  of  a  wln^  with 
a  secondary  control  surface  It  is  necessary  to  replace  a  by  (a  +  6) 
and  <5  by  —6  In  the  formulas  found  for  a  wing  with  a  flap.  Carrying  out 
this  operation  in  Formula  (3.16),  we  will  obtain 

7  = /»?  («  +  fl)  (g  »B  T -r ''wfl  tc  e  (2fcA -f  6>  tg  T).  (3.24) 


Hence  we  find  that 


Z  —  Ic  r.S^  -i-  tc  e5„, 


(3.25) 


Op  * 

where  S  =  [B  tan  y  -(2bh  +  b  tan  y) ]  is  the  area  of  the  secondary 
control  surface. 

For  the  coefficient  of  force  (3-25)  we  will  obtain  the  following 
expression; 

c,  VM*—i 


4a 


-(‘■■-•ft)- 


(3.26) 


Using  Formula  (3.18),  we  will  find 

M-=4galgea*<gTY- ^  + 

-  (2bh  -i-  fcMg  T)  Y -J 

For  the  coefficient  of  the  moment  (3.27)  we  will  obtain  the  fol¬ 
lowing  expression: 

■+T[TT-r(-j)x+TT^]-  (3-28) 


.9  - _ 


4a  3 

The  formula  for  the  determination  of  x^^  ^  takes  the  following 


form; 


-  40  - 


PTD-TT-62- 1694/1+2 


where 


^x, 


HA  = 


.  ^ 

1 

5,.;, 

±T, 

•  *« 

3  i 

B  ) 

6 

5k 

J 

(3.30) 


Carrying  out  the  substitution  In  (3.22),  we  will  obtain 

T  T  [ -ifr + 7?^ 

For  the  coefficient  of  the  moment  (57)  we  will  obtain  the  follow 


tg«e  tgr]  (331) 


3. 32) 


Ing  expresslonsjvV^^  ,  j,,  2ri«e  _ tgr  »geT 

\  a  /  3  jf  [tgT  ^  —  tgr  j 

Dividing  (3-32)  by  (3.26),  we  will  find  the  dimensionless  coord¬ 
inate  of  the  center  of  pressure  of  the  forces  acting  on  the  half  wing 
about  the  half -span. 

According  to  the  found  formulas  for  c^  and  ^  we  carried  but 
the  calculations  which  were  then  compared  against  the  experimental  re 
suits.  The  comparisons  showed  that  the  calculated  data  are  In  good 
agreement  with  the  experimental  within  the  following  range  of  numbers 
2  <  M  <  4. 1. 
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[List  of  Transliterated  Symbols] 

U.fl  a  ts.d  =  tsentr  davlenlya  =  center  of  pressure 
3  =  e  =  eleron  =  aileron 
K  =  k  =  krylo  =  wing 


-  41  - 


PTI>-  TT-  62- 1694/1+2 


39  3  =  z  *s  zakrylok  =  flap 

39  3.n  =  z.p  =  zakrylok  predkrylok  =  flap  and  other  control 

surface 
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